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Abstract 

We consider various forms of the mass term that can be used in 
the Skyrme model and their implications on the properties of baryonic 
states. We show that, with an appropriate choice for the mass term, 
without changing the asymptotic behaviour of the profile functions at 
large r, we can considerably reduce or increase the mass term’s contri¬ 
bution to the classical mass of the solitons. We find that multibaryon 
configurations can be classically bound at large baryon numbers for 
some choices of this mass term. 


1 Introduction 

The Skyrme model has enjoyed a lot of interest ever since it was realised that, 
although it is a nonlinear theory of pions, it is also an effective theory of low 
energy nucleon interactions. In fact, it may also provide a new approach to 
nuclear physics; as the lowest states of the model, corresponding to higher 
baryon numbers, are expected to provide a classical description of nuclei. In 

*email: kopelio@al20.inr.troitsk.ru 
lemail: b.m.a.g.piette@durham.ac.uk 
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the Skyrme model approach the baryon number is identihed with the soliton 
number. 

Multiskyrmions are the stationary points of the static energy functional 
which, in natural units of the model, is given by 


E = 


127r2 



diUU-^Y + ^Tt 

1 2 

J 16 L 


diUU 


-1 


djUU-^ 


1 2 


d^x. ( 1 ) 


where U{x) G SU{2) and x is in units of 2/{F.„e). 

Most of the phenomenological applications of the Skyrme model, espe¬ 
cially to the study of the nucleon or hyperon properties, included also the 
pion (or kaon, or D-meson) mass term in the Lagrangian chosen in the sim¬ 
plest possible form (see eg Adkins and Nappi ^J). In particular, the kaon 
mass term has to be added to describe the mass splittings within the SU (3) 
multiplets of baryons: octet, decuplet, antidecuplet, etc |2]. However, the 
role of the mass terms in multiskyrmion conhgurations, especially at large 
baryon numbers, has not been investigated in much detail; the theoretical 
work performed so far has involved mostly the Skyrme model in which pions 
are massles [i.e. given by the Lagrangian above). It is only very recently 
that some attention has been paid also to the effects associated with the pion 
mass for large B conhgurations inmi; one of the effects being the exponential 
localization of multiskyrmions. In particular, it was stressed that the contri¬ 
bution of the mass term can change the binding properties of large B classical 
conhgurations and, in particullar, their decay properties into conhgurations 
with smaller H-numbers jl]. 

In most of these approaches the pion mass term has been introduced via 
the addition to (CD of the following term 

L ( 2 ) 

where m is related to physical pion mass fi-jr ~ 138 Mev by the relation 
m = 2/i7r/(i^Trc), where ~ 186 Men is the pion decay constant, taken usu¬ 
ally from the experiment, and e is a Skyrme constant*. The appearance of 
the mass term in ehective held theories was discussed, e.g. in jS]. Although, 

*In ^ the masses of the nucleon and the A(1232) isobar were fitted using an SU{2) 
quantization procedure and, as a result, the authors obtained = 108 Men, e = 4.84, 
but these values did not allow to describe the mass splittings within SU{3) multiplets of 
baryons. The approach of ^ has been revised and another set of parameters is widely 
accepted now. The baryon mass splittings are described with the experimental values of 
Ftt, Fk and e ~ 4.1. For these parameters the absolute values of the baryon masses are 
not fitted because they are controlled by the loop corrections, or the so-called Casimir 
energy which, for the baryon number 1, was estimated in [S]. 
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the effects associated with the pion mass are small for the small values of this 
mass, they increase if either the baryon number or the pion mass are larger. 
For massless pions all the known minimal energy multiskyrmion configura¬ 
tions have a shell-like structures. These field configurations were obtained 
in both numerical simulations and in studies involving the so-called ‘rational 
map ansatz’. In the rational map ansatz one approximates the full mul¬ 
tiskyrmion field by assuming that its angular dependence is approximately 
described by a rational map between Riemann spheres. This approximation 
was, hrst of all, shown to be very good in a theory with massless pions and 
it was later extended also to massive pions - where the agreement was again 
shown to be very good. 

Given these observations it is extremely important to have the right (cor¬ 
rect) mass term. The problem, however, is that the mass term is very non¬ 
unique and the expression m is only one of many that can be used. Indeed, 
the origin of the chiral symmetry conserving and chiral symmetry breaking, 
or the mass terms considered in [^, may have a very different nature. So we 
have decided to reexamine this issue further and to look at mass terms other 
than 0 and see what effects they have on the properties of multiskyrmion 
conhgurations. 

In the next section we discuss various choices of the mass term, pointing 
out what is hxed and what can be changed, and in the following sections 
we look at some simple examples of such mass terms. Expressions for the 
static energy of skyrmions and some dehnitions necessary for the description 
of multiskyrmions within the rational map approximation are presented in 
sections 3 and 4. Our numerical results are presented in section 5 and the 
analytical discussion useful to establish asymptotic behaviour is presented in 
section 6. We finish with a short section discussing our conclusions and ideas 
for further work. 


2 Mass Terms 

To consider the mass term we hrst note that the pion helds if = (tti, 712 , 713 ) 
are given by f/ = a + iw ■ f, where r denotes the triplet of Pauli matrices and 
a is determined by the constraint cr^ -|- vf • vf = 1. 

Then the square of the pion mass is the coefficient of the expansion of the 
mass term in powers of vf ■ vf; in fact it is the coefficient of the lowest term 
i.e. vf ■ vf in this expansion. In the case above we have 

m^Tr(l — U) = 2(1 — a) ~ rn? vf • vf -|- ..., (3) 

where +... stands for further powers of vf ■ vf to be interpreted as pion inter- 


3 


action terms. So the mass of the pion held is proportional to m, since the 
canonical mass term in the lagrangian is —■ 7f/2. 

However, 0 is not the only term we can use as the pion mass term. It 
is clear that we can multiply (1 — t/) in Q by any function of U which in 
the limit > 1 reduces to 1. Thus we could multiply it by, say, 

In fact, a little thought shows that, instead of U in 0 , we can take 

/ OO 

g{p)UPdp (4) 

-OO 


where 


dip) dp = 1 and 


g{p)p^dp = 1. 


(5) 


J—OO j—OO 

The usual choice then corresponds to 

g{p) = dip-l). (6) 

As the second condition in © can be eliminated by redehning the coef- 
hcient m? in m we see that a more general mass term is given by 


1 

127r2 


Am^ 



/ OO 

gip)U^dp 

-OO 


d^x, 


(7) 


where _ 

/ OO 

g{p)p^ dp. 

-OO 


( 8 ) 


3 B=1 Skyrmion 


Consider hrst the case of one Skyrmion. The single Skyrmion has the hedge¬ 
hog form 

U = exp(i/(r)r • r), (9) 

where r is the unit vector in the f direction and /(r) is the radial prohle 
function which is required to satisfy the boundary conditions /(O) = tt and 
/(oo) = 0. 

Putting 0 into the energy functional we hnd that the energy of the held 
is given by 


E 


1 


Svr Jo 
2Am?r‘^ 


+ 2(/2 + l)sinV + 

roo 

1 — / gik) cos{kf ) dk 

J —OO 



( 10 ) 


where A is given by 0. 
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Thus, for the minimal field, f{r) satisfies the equation 


/[r^ + 2sin^/] + 2r/ + 2(/^ — 1) sin / cos / 

2 2 I-oc 9{k) sm{kf)kdk 

— 777 T - 

IZo9ik)k^dk 


2 sin^ / cos / 


0 . ( 11 ) 


We have investigated several classes of such functions: 

• g{k) = 5{k — p) for several values of p. 

The cases of even or odd integer values for p have also been investigated 
analytically. Moreover, it is easy to notice that the mass term for p > 1 
is smaller than that for p = 1, since (1 — cos(p/))/p^ = 2sin^(p//2)/p^ 
and sin^(p//2)/p^ < sin^(//2) for p > 1. 


• g given by a Gaussian centered around p = 1, or around p = 0. 
In the latter two cases we have taken 


dip) = ^exp (-A(p - 1)2) and g{p) = ^exp (-Ap^) . 


The corresponding expressions for the energy are given by 

E = 7rl>'V= + 2(h+l)»i»V + ^ 


and 


Sir Jo 


2m^r‘^ 


^ + 1 + ^ 
v/;^A ^ ^ ^ 2A L 


f2 

1 -cos/exp(-^) 


>dr, 




Stt Jo 

r p • 

+ Am^r^X 1 — exp(——) 

4A 


>dr,, 


respectively. 


( 12 ) 


( 13 ) 


( 14 ) 


4 Multiskyrmions 

For multiskyrmion fields we use the rational map ansatz of Houghton et ah 
[2j. The ansatz involves the introduction of the spherical coordinates in R^, 
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so that a point x G is given by a pair (r, where r = |x| is the distance 
from the origin, and is a Riemann sphere coordinate giving the point on 
the nnit two-sphere which intersects the half-line throngh the origin and the 
point X, i.e., ^ = tan(|)e®‘^, where 9 and (p are the usnal spherical coordinates 
on the nnit sphere. 

Then one observesjS] that a general SU{2) matrix, U, can always be 
written in the form 

U = exp{if{2P - /)) (15) 

where / is real and P is a 2 x 2 hermitian projector i.e., P = = Ph The 

rational map ansatz assnmes that the Skyrme field has the above form and, 
in addition, that / depends only on the radial coordinate, i.e., f = /(r), and 
that the projector depends only on the angnlar coordinates,he., P(^,^). 

The projector is then taken in the form 


P 


f 

IfP 


( 16 ) 


where f(.^) is a 2-component vector, each entry of which is a degree k poly¬ 
nomial in Incidentally, given the projective natnre of f one can also nse 
the parametrization f = (1, P)*, where P(0 is the ratio of P = ^. 

For P = 1 this ansatz reprodnces the one Skyrmion field confignration 
discussed in the last section, while for P > 1 the ansatz (USD is not compat¬ 
ible with the equations which come from o, so the ansatz cannot produce 
any exact multi-skyrmion configurations. However, as was shown in many 
papers, see eg: 0, i, it gives approximate field configurations which turn 
out to be very close to the numerically computed minimal energy states. To 
do this one selects a specihc map f and puts it into the Skyrme energy func¬ 
tional (P). Performing the integration over the angular coordinates results 
in a one-dimensional energy functional for /(r) which has then to be solved 
numerically. 

Hence, if we write f = (1, P)*, then the Skyrme energy is 


E 


- rV'^ + 2P(r + l)sinV-f J 


-sinV 


Stt 
2Aw?r‘^ 


dr 


1 - 


g{k) cos{kf) dk 


dr, 


where X denotes the integral 


J = 


1 

dvr 


/1 + iep 

dR 

y 2 i d^di 

U + |PP 

di 

) (1 + ieiT 


(17) 


( 18 ) 
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The values of X have already been calculated; in what follows we take our 
values from [H] . The equation for the prohle function / is very similar to the 
equation HU of the last section except that the coefficients of terms involving 
sin / cos / are multiplied by B in the first term and X in the second. 


5 Numerical results 


We have looked at the values of various quantities for different choices of p 
(taking g{k) = 5{k — p)), and also at some Gaussians. The Gaussian cases 
were not particularly illuminating so here we discuss only the cases of fixed 
values of p. 

Note that when p —0 we have a nontrivial contribution of the mass term. 
This involves taking the limit p —> 0 of the expression in and then its 
last line becomes We can also consider the limit when m —> 0 which 

corresponds to the massless Skyrme model. 

We present our numerical results in hgures ITIIHI in which we plot the 
normalized energy 


En 


E(B) 

BB(1) 


( 19 ) 


and the shell radius, as a function of B, for several values of the mass 1 and 
for p from 0 to 5. 

The normalized energy (m is a dimensionless energy which describes the 
binding of the configuration by comparing it to that of the B = 1 solution. 
Note that when En > 1 the B multiskyrmion configuration has an energy 
larger than the energy of B single skyrmions thus showing that this config¬ 
uration is unstable. The jagged curve near the origin is caused by the value 
of X which varies a lot when B is small. When B > 22, we have taken 
X = 1.285^ (see jHI) and so the curves are smooth. 

We also present the radius of the solutions dehned as follows: 


R 


J rE{r)r^dr 
f E{r)r^dr ’ 


( 20 ) 


where E{r) is the radial energy density. 

Looking at the figures we see that, as B gets very large, the normalized 
energy converges to a finite value when p is even, but that it slowly diverges 
when p is odd. This was also observed by Battye and Sutcliffe jlj in the case 

Tn our numerical calculations we take for pions m = = 0.36192, for kaons rUst = 

1.29996 and for the charmed mass scale rrich = 4.130964. In what follows, in the text and 
in the captions, we refer to those values as m = 0.362, 1.300 and 4.131 respectively . 
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of p = 1. When m = 4.131 (Fig. [T)) the curve for p = 3 crosses the value 
of En = 1 at i? ~ 365. For smaller value of m, the curves cross this bound 
state threshold for much larger values of B. This is not surprising, as for odd 
p, the mass term adds a non vanishing contribution to the energy density 
from the region where / = tt and so to make a large shell conhguration one 
needs to put an increasing amount of energy inside the shell. When B is 
large, this becomes energetically too expensive and the conhguration is not 
a bound state anymore. 

We also see that, for a given parity of p, the energy at a given value of B 
decreases as we increase p by a multiple of 2, i.e. Enp{B) > Enp+ 2 {B)). 

The plots of the radius also show that when p is odd the shell is smaller, 
but otherwise, for a given parity of p, the radius increases with p. On the 
other hand, for hxed values of p and B, the radius decreases when the mass 
m increases. This is exactly what one expects for odd values of p as the 
energy density inside the shell is non zero, but it is also true for even p, i. e. 
the radius of the shell decreases with the increase of the mass. 



0.8 


100 200 300 400 500 


100 200 


300 400 500 


B 
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a 
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Figure 1: Normalized energy (Uni) (a) and radius (pr)|) (b) of multiskyrmion 
conhgurations for m = 4.131, E{1) = 2.056 


One property worth investigating for these low energy conhgurations is 
their ability to decay into two or more shells of smaller baryon charges. To 








a b 

Figure 2: Normalized energy (0 (a) and radius (b) of multiskyrmion 
configurations for m = 1.300, E{1) = 1.486. 

do this we have computed the derivative of the energy with respect to B and 
compared thus obtained values with the energy per baryon of some small B 
conhgurations of low energy (typically B = 2, B = 4:, B = 7 and B = 17). 
When the value of the derivative is larger than the energy per baryon of 
other configurations, it implies that the larger conhguration can decay into 
two shells. In Fig. I^we see that, when m = 0.362 and p = 1, all large shells 
{B > 50) can decay into shells with B = 17, B = 7 and B = 4 but can only 
decay into a B = 2 ii B > 250 and i? = 1 if i? > 380. When p is even but 
non zero, the normalized energy decreases as B increases implying that the 
binding energy of the conhgurations increases with B. This in turn implies 
that the conhgurations never decay into smaller shells. 

We summarize our observation in Tabled where we have given the thresh¬ 
old value corresponding to several decay modes. The threshold values are the 
values of B above which the decay is always possible, but sometimes some 
conhgurations with lower values of B (smaller than 32) can decay in such a 
mode too (but only for very special values oi B). The values given with a 
‘>‘ refer to the values obtained by comparing the energy of the conhguration 
directly (low B) instead of using the derivative of the energy. 

In Fig. El we present the normalized energy (unD as a function of m for 
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Table 1: Decay of low energy configurations into sub shells. Each column 
corresponds to a decay mode and the baryon charge given corresponds to the 
threshold from which the decay is always possible. 
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a b 

Figure 3: Normalized energy (0 (a) and radius (b) of multiskyrmion 
configurations for m = 0.362, E{1) = 1.274. 


various values of B. We see that for p = 0 and p = 1 the energy increases 
rapidly with m and that very quickly the conhgurations become unstable. 
When p > 1, the normalized energy decreases for small m when m increases, 
then it reaches a minimum and hnally it increases with m. The value of 
the mass for which the solution is the most bound depends on p and on the 
baryon number. 

In the next section we present an analytical description of the multi¬ 
skyrmion conhgurations which will explain some of the features we have 
observed numerically. 


6 Approximate analytical treatment 

It was shown in [Hj that many properties of multiskyrmions, including their 
classical mass, spatial distribution, moments of inertia, etc., can be described 
with a good accuracy using a relatively simple power-step (or “inclined” step) 
approximation of the prohle function. A similar approach was also used 
successfully to describe “baby”-skyrmions in the 2-1-1 dimensional Skyrme 
model [na [mini CHI. This approximation also turns out to be useful for the 
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Figure 4: Derivarive of the normalized Energy with respect to B and the 
normalized energy oi B = 2, B = A, B = 7 and B = 17 for m = 0.362 and 

p = 1. 

study of the asymptotics of our massive multiskyrmion held conhgurations 
for large values of the baryon number B. 

Let us consider hrst the large r asymptotics of the prohle function. Clearly, 
this asymptotic behaviour is governed by the 2nd derivative term and the 
mass term in the Lagrangian. The Euler-Lagrange equation then becomes 
asymptotically 

(rY)'=(2i? + mV)/, (21) 

and so, if 3> 2B, we have 2rf' + r^/" = which has the asymp¬ 

totics / ~ exp{—mr). 

For the values of B in the region of r where < 2B the prohle / 

has a power behaviour, and it is in this region that most of the mass and 
most of the baryon density of the multiskyrmion is concentrated |21, while 
the exponential tail of the prohle function gives only a small correction to all 
such quantities and so can be neglected. Thus if we can neglect the ~ m? 
term on the right side of we obtain the power law f ^ r , As 
we shall see, the dominant range of r is always such that we can make this 
approximation, at least for pions and kaons. 
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Denoting 0 = cos / and taking for g{k) = 6{k — p) where p is an integer, 
the energy of multiskyrmion can be written as 


M 


Stt J I (1 — 0^) 
250'^ +1 


[rV'^ + 2S(1 - 02)2] 

+ p{(j))P \dr, 


( 1 - 02)2 


( 22 ) 


with 0 changing from —latr = Otolatr—>oo. The hrst part of is the 
second order term contribution while the second term is due to the Skyrme 
term. Note that at hxed r = vq the 4-th order term is exactly proportional 
to a 1 -dimensional domain wall energy widely discussed in the literature, see 
e.g. ^3]. The function 4/^(0) = (1 — cos(p/))/p2 can be written explicitly for 
each p\ Ti = 1 - 0, T 3 = (1 - 0)(1 20)2/9 < Ti, T 2 = (1 - 0^)/2 < Ti, 

4/4 = 02(1 — 02)/2 < 4 / 2 , etc. Also it can be shown that 4/3 < 4/2, tt '4 < 4 / 3 , 
and it follows immediately, for any p, that 4 / 4 ^ < 4 / 3 ^ < 4/2p < 4/p, etc. The 
functions 4/p and the whole mass term have different properties for odd and 
even p and so, for this reason, these two cases will be considered separately. 

It is possible to rewrite the second order term contribution in ()22j) as: 


M( 2 ) = 


StT y 1(1 — 02) 


-|2 >1 

0'- \/^(l - 02)/r +2rv^0'Wr, (23) 


and similarly for the 4-th order Skyrme term. Next we observe that if 0 
satishes 0' = \/2B{l — 4>^)/r, a large part of the integrand in M*^ 2 ) vanishes. 
Therefore, it is natural to consider a function 0 which satishes the following 
differential equation jS]: 

^' = 4(1-^^), (24) 

where 6 is a constant. A solution of this equation, which satishes the bound¬ 
ary conditions 0 ( 0 ) = —1 and 0 (cx)) = 1 , is given by: 


0 (r, ro, h) 


(p/Po)b _ 1 

{r/roY + 1 


(25) 


where tq is the distance from the origin to the point where 0 = 0 and at 
which the prohle f = 7i/2. tq can be considered as the radius of the multi¬ 
skyrmion. Both b and ro are arbitrary at this stage; they will be determined 
later by means of the mass minimization procedure. Note that the radii of 
distributions of baryon number and of the mass of the multiskyrmion are 
close to Tq. Let us point out that our parametrization (123) is very accurate 
as, in the Skyrme model with the usual mass term, as shown in j2j, the 
masses and other characteristics of multiskyrmions are described by such a 
parametrization to within a few %. 
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6.1 Odd powers, ]? = 1, 3,... 

Consider first the case of p = 1. Then, nsing 03 and we find that the 
soliton mass is given by 




1 

Stt 



+ 2s)(1-02) 


+ 

+ 


2r^m^(l — 4>)\dr, 


2 


+ 


(26) 


where 0 is given by (j2Sl) and where we shonld take m = 0.362 for the pion 
case and m = 1.30 for kaons, etc. 

Given the form of 0 the integration over r can now be performed nsing 
the well known expressions for the Euler-type integrals, e.g. 


dr 


Tiro 


lo l-|-(r/ro)^ 6sin(7r/6)’ 


if 6 > 1, and, more generally jS], 


( 27 ) 


(r/roYdr ^ Trrp 

Jo P+ir/roY 6 sin[7 r(lc)/ 6] ’ 

with /9>0, 6 >l-t-c, c>—1. Differentiating with respect to [3 allows us 
to get the integrals with any power of 1 -|- {r/ tq)^ in the denominator. Thus 
we can derive the following expressions for the integrals of 0 given by 


y (1 — 0^) dr = 


dvrrn 


2\2 


(1 - 0^) 


dr = 


87r(l - 1/&2) 


62sin(7r/6)’ J 3ro6^ sin(7r/6) ’ 

and other examples useful for the calculation of the mass term. 


( 29 ) 


fd,M-d>Vir= ( 30 ) 

J bsin[6TT/b) j b^sm[6n/b) 


{Sn/b) 


The expressions fl2hl3()j) allow us to obtain the mass of the multiskyrmion 
held in a simple analytical form as a function of the parameters b and ro (in 
units 37r^F,r/e): 


M{B, ro, b) = a{B, b)ro + (5{B, b)/ro + 6{b) r^ (31) 


where 


a = (6^-I-8i?)/(36^sin(7r/6)), 

0 = 4{Bb^+ 21){l-l/b^)/{9b^sm{7^/b)), 

6 = 4m^/(36 sin(37r/6)). (32) 
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The mass term contribution is proportional to the volume of the multi- 
skyrmion, ~ Tq, as expected on general grounds, multiplied by the corre¬ 
sponding flavour content of the skyrmion h Next we minimize ea with 
respect to tq and obtain, in a simple form, the precise minimal value of the 
mass 

n^min _ 

M{B,b) = ^-{^a^ + 12Sl3 + 2a) (33) 

o 

where the value of tq is given by: 


min 

^0 


{B,b) 


Va^ +12/3(5 
6(5 



( 34 ) 


Eq. (j31l give the upper bounds for the mass of the multiskyrmion state, 
because they are calculated for the prohle (12511 which is different from the true 
prohle to be obtained by the true minimization of the energy functional (uni) 
with the mass term included. At large values of B the power b is also large, 
b ~ ^/B, as we shall see, and o; ~ (6 -|- 8i?/6)/(37r), (3 ~ A{Bb + 2X/6)/(97r), 
6 ~ 4m^/(Ovr). 

The structure of m remains the same for values of p different from 
1, except for the case of even p which will be considered separately. For 
p = 3, 5, etc. one must perform the substitution 6 —>■ 5/p^, and the volume 
contribution is reduced by a factor 1/p^. The energy (|HH|l can be simplihed 
and analyzed in two different cases, small m or 6, when 12(5/5 -C (which 
we will call in what follows the small mass approximation, or SMA), and 
in the case of large m or large B when 12(5/5 /:§> o;^, which we will call the 
large mass approximation, or LMA. Note, that at large i5-numbers, a r\j y/B 
and (3 ~ B\fB] therefore, when B is large enough, the latter inequality can 
always be satisfied: it reads then, approximately, m?y/B S> 1. 

Let us consider first the latter case of large (35 (the LMA case). Now we 
can neglect the term ~ o? in the square root of ()33l34j) . and obtain 


M{B,b)^^^[(3H 


1/4 


3 a/3V/2 

^ 4(/5(5) . 


(35) 


and 


. 1,1.. L. 

I n — 


£xl/4 

36 ) 


-I 


1 y/2 


(36) 


Mhere is a difference of principle between the pion and kaon (or other flavoured meson) 
fields included into the Lagrangian. The mass of flavoured mesons enters into the classical 
soliton mass multiplied by a corresponding flavour content which is always smaller than 
1 and even smaller than 0.5 for a rigid or soft rotator quantization scheme. So, when we 
take the masses m ~ 1.30 for the strangeness or m ~ 4.13 for the charm we establish the 
scale of the mass for these flavours, not more than that. 
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B 

p=l (num.) 

p=l (SMA) 

p=3 (num.) 

p=3 (SMA) 

1 

1.2740 

— 

1.2576 

— 

40 

1.1519 

1.1497 

1.0990 

1.0935 

100 

1.1734 

1.1760 

1.0991 

1.0982 

200 

1.1973 

1.1956 

1.1023 

1.1034 

300 

1.2144 

1.2037 

1.1053 

1.1073 

400 

1.2279 

1.2057 

1.1080 

1.1106 

500 

1.2392 

1.2038 

1.1104 

1.1133 


Table 2: Energy per baryon for m = m^r = 0.362. The analytical calculations 
are made in the small mass approximation (SMA) according to (HTl) . 


It is clear that the minimum value of the mass is reached at the minimum 
of (3 {6 does not depend on b when b is large, and the correction term in the 
square bracket has little influence on the position of the minimum), which is 
equal to = 8-\/2i?J/(O tt) at b = Then 


M{B) ~ 




97rV3y 


m^/\2Bl) 


3/8 


1 + 


3^3 (J +4^2) 
^y/2m{2IBfl^ 


(37) 


Since X 52 (strictly, X > 52 [Zj), we establish the following scaling law at 
large B-. M{B) ~ r{B) ~ Numerically, we have for 

p = 1 (X = 1.28i?2 in these estimates) 


M, , ^ UR f 1.1982\ 

(p = 1) ~ 0.59395v^5^/® I 1 + —^ 1 


B 


(38) 


For other odd p, dividing the volume contribution to the mass term by p2, 
we obtain 


M, , 0.59395 1.1982 \ , , 

-(p).—(39) 

The absolute lower bound for the energy which follows from dSi obviously 
does not depend on p. These estimates can be improved further: the 0{a^) 
terms in the expansion of the square root in dSl can be included; the surface 
contributions to the mass term, besides the volume-like one, can be calculated 
(this may be important for higher p since the volume contribution decreases 
like ~ 1/p^); the shift in the position of 6™" could also be taken into account. 

As the ratio M{B)/B ~ we conclude that, for large B, the shell 
conhgurations will not form a bound state. This conhrms what we have 
observed numerically. The second order term of the initial Lagrangian makes 
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B 

p=l (num.) 

p=l (LMA) 

p=3 (num.) 

p=3 (LMA) 

p=3 (SMA) 

1 

1.4860 

— 

1.3842 

— 

— 

40 

1.4525 

1.4675 

1.1893 

1.3018 

1.1701 

100 

1.5375 

1.5553 

1.2097 

1.3032 

1.1968 

200 

1.6181 

1.6351 

1.2360 

1.3157 

1.2056 

300 

1.6714 

1.6875 

1.2554 

1.3276 

1.1982 

400 

1.7120 

1.7273 

1.2710 

1.3381 

1.1820 

500 

1.7450 

1.7596 

1.2841 

1.3474 

1.1603 


Table 3: Energy per baryon for m = = 1.300. The analytical estimates 

are made according to (IHH]) (LMA) and, for p = 3, also in SMA. 


B 

p=l (num.) 

p=l (LMA) 

p=3 (num.) 

p=3 (LMA) 

1 

2.0558 

— 

1.7370 

— 

40 

2.1778 

2.1352 

1.5160 

1.4877 

100 

2.3619 

2.3436 

1.5839 

1.5805 

200 

2.5304 

2.5216 

1.6589 

1.6643 

300 

2.6398 

2.6344 

1.7111 

1.7191 

400 

2.7222 

2.7185 

1.7516 

1.7607 

500 

2.7888 

2.7861 

1.7850 

1.7945 


Table 4: Energy per baryon for m = rrich = 4.131. The analytical estimates 
are made in the large mass approximation, ()38I39|1 . 


a small contribution in the large mass regime; thus the Skyrme term and the 
mass term approximately balance each other, and the mass term gives 1/4 
of the total mass, by the Derrick theorem. The difference between the cases 
p = 3, 5,... etc. and p = 1 resides in the fact that for larger p this “large 
mass term regime” is reached at higher values of the baryon number. One 
of the properties of multiskyrmions in this regime is that the average energy 
density does not depend on B, pM ~ , or, in ordinary units, pm ~ plFn 

(which does not depend on the Skyrme parameter e). The energy density 
in the shell can be estimated as well; we get pshell ~ which grows 

when the baryon number increases. And, as it has been previously discussed 
in the literature, the transition to other types of classical conhgurations, like 
the skyrmion crystals, may become possible at high values of B. 

When the mass m is small enough, as for the pion, the expansion in 
12j35/a^ can be made, and one obtains the reduction of the multiskyrmion 
size Tq: 


ro 






-jV4 




(40) 
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and the increase of its mass 


6M — Mm=o 


(35 


9(35 






1 - 


(41) 


r]B = Vl/{2B + VX), M^=o ^ AB^{2 + eB)/(37r) 0, and at large B, 

Mm=o/B ~ 1.0851, rjB — 0.3613 if we take ^b = — 1-13137 - constant 

valne, according to p. Here, we have nsed also the observation that, at large 
5, = 2J^/^ P and 

As expected, the size of the mnltiskyrmion state decreases with increasing m 
while its mass increases, and these changes become very large for very large 
B and/or m. 

For p different from 1 the snbstitntion —>■ (p^ should be made in 

mm and the following relation can then be obtained for any pair of odd 
p’s. Pi and p 2 . 


Mb{Pi) - Mb{p2) ^ ^ 71 _ 1 Wi/4„^ .42) 

Mb 9 Vp? pV ^ ’ 

Numerically this works well for p = 1 and p = 3, see Table l2im for larger p's 
the agreement is less good but then, apparently, other contributions to the 
mass term, besides the volume-like one, should be included. 

In Tables 121111 we present several values of the energy per baryon obtained 
from (Table P and and (P?|l (Tables mH), and compare them with 
the values obtained numerically. We see that our analytical approximation 
works very well when the mass is small and H-numbers not too large (pions 
case. Table ED, or when it is large, as for the charm, but it does not work 
so well for intermediate values of the mass. It also works better for p = 1 
than for p = 3. The case p = 3, presented in Table P is of special interest: 
the LMA improves when increasing the baryon number but is still not as 
good as for p = 1, whereas SMA becomes worst when B increases and also 
is not perfect at small values of H. To improve it, the following terms in the 
expansion (ED should be included. 

For not very large values of m the structure of the mnltiskyrmion at large 
B remains the same: it is given by the chiral symmetry broken phase inside 
a spherical wall where (on this spherical shell) the main contribution to the 
mass and topological charge is concentrated PEI. The value of the mass 
density inside this wall is dehned completely by the mass term with 1 — 0 = 2 
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and decreases with increasing B while the mass density of the shell itself is 
constant [3]. The baryon nnmber density distribntion is qnite similar, the 
only difference being that inside the spherical wall it vanishes. 

If, for some physical reasons, we wonld nse a distribntion over p in the 
Lagrangian, as discnssed in Sections 2,3, the analytical expression for the 
mnltiskyrmion energy can be obtained qnite analogonsly. Let ns pnt p = 
Po + Ap, where po = 1, 3, etc., and Ap is assnmed to be small. Then, taking 
into acconnt the changes in the volnme contribntion to the mass term, we 
get instead of (inni): 


M 

'b 


(P) 




0.59395 

VPo 




1 


X 


1 +Po 


Ap (Ap)2/ 2 

— - — -I TT 

2po 16 V 


Po 



\ 1.1982 

JmB^B ’ 


( 43 ) 


and the averaging over any distribntion g{p), as snggested by ()4l5ll2j] can be 
easily performed. 

It is also possible to consider, in a similar way, the case of small valnes 
of p, near p = 0. In this case we have (1 — cos(p/))/p^ ~ f^{l — p^/^/12)/2, 
and / ~ TT inside the mnltiskyrmions. Evalnations similar to those at the 
beginning of this section show that the energy per nnit i?-nnmber, in the 
large mass regime, is given by: 

— ~ 0.74441-\/mi?^'^® f 1 — 

B V 




2^2 


1 + 


0.7628 / ^ p"7i 


mB^B \ 


1 + 


24 


(44) 


Obvionsly, at large enongh valne of the mass, this is somewhat greater than 
the energy given by ()38|1 . The expression il can be compared with onr 
nnmerical resnlts for p = 0 presented in Fig. [T] and Fig. |21 We note the 
agreement to within an accnracy of abont (5 — 8)% for the largest valnes of 
m and B. The integration over any distribntion in p, g{p), near p = 0, as 
presented in (M, can also be easily made. 


6.2 Even power p = 2, 4,... 

In the case of even p he., p = 2, 4... the volnme contribntion to the energy 
density is rednced becanse 1 — cos(p/) ~ 0 inside the mnltiskyrmion where 
the profile / — vr. Dne to the dependence of the mass term on the parameter 
b and dne to the connection between tq and b this case is very different from 
the case of p = 1, 3.... However, we can still write 


M{B,ro,b) 


~ a{B, b)ro + 


P{B,b) 

ro 


b 


(45) 
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This expression coincides with (EH), except for the mass term where an ad¬ 
ditional factor 1/6 appears and 6' is different from 6. For p = 2 we have 
6'{p = 2) = Am?/{h sin(37r/6)) ~ 4m^/(37r) at large 6 while for p = 4 an ad¬ 
ditional small factor appears as 5'{p = 4) = 4m^(l -|- 18/6^)/(36 sin(37r/6)) ~ 
4 m^/(97r). It is not easy to hnd the general expression for larger p; i.e., the 
expression for 5' which decreases with increasing p (see the discussion after 

m)- 

In general, we proceed as for odd values of p and after minimizing with 
respect to ro we obtain 


^j^rrnn _ 

M(B, 6) = —+ i2S'p/b + 2a) 


and the value of ro 


r7"(5,6) = 


\JoP‘ + 125'(3/b — a 


66'/b 


1/2 


(46) 


(47) 


The main difference from the previous case is that, at large values of -B, 
the quantities ~ -B and 6'13/b ~ B, i.e., they are of the same order of 
magnitude, since 6 ~ \fB. At large enough 6 or i? we have 12(56'/{ba^) ~ 
0.35 for pions, 4.5 for kaons and ~ 46 for charm. 

Let us discuss hrst the large mass case when we can take 


^0= + 121';3/6 ~ 2,j3S'p/b 


(48) 


and 


4 


. a ( 3b \ 


(49) 


The minimum is reached at 6 ~ 2^X/B and we have, recalling that at large 
B and p = 2, A ~ 4m^/(37r), 


M{B,p 


2 ) ~ 


3v/3 7r6V4 


\^3V3iiB±2l^ 

4:\/2m 


(50) 


at ro ~ yi^/m. Note that and, at large B, it is constant within 

the rational map approximation. 

Numerically, for p = 2, we obtain from 


M 

~B 


(p = 2) ~ 0.66612 \/rn 



0.8877 

m 


(51) 
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m 

p=2 (num. B=500) 

p=2 (approx) 

p=4 (num. B=500) 

p=4 (approx) 

0.362 

1.0879 

1.1016 (SMA) 

1.0881 

1.0908 (SMA) 

1.300 

1.2186 

1.2048 (SMA) 
1.2781 (LMA) 

1.1362 

1.1475 (SMA) 
1.2596 (LMA) 

4.131 

1.6236 

1.6448 (LMA) 

1.3519 

1.4116 (LMA) 


Table 5: Asymptotic values of the energy per baryon for p = 2 and p = 4. 
The analytical calculations correspond to dnH), SMA, and (I51I52I) in LMA. 


and for p = 4 

M, , 1.5375\ , , 

— p = 4 - 0.50614 l +- . 52 

B \ m J 

In Table 0 we present a few values of the asymptotic energy obtained 
from (IMD, SMA, and from m and (IS2D in LMA and compare them with 
the values obtained numerically for B = 500. In our calculations, for large B, 
we have again used the value = 1.13137, P, and 6' ~ 4m^/(O tt) for p = 4. 
There is a good agreement between our numerical results and our analytical 
approximation values when the mass is small (pions), or large (charm scale), 
and not so good for the intermediate value m = rUst where we give both 
the SMA and LMA results. Our approximation also works better for p = 2 
than for p = 4. Nevertheless, analytical approximations work better for odd 
p at large m and B, than for the even ones. It is possible to improve the 
analytical estimates although, for even p, the estimate of the preasymptotic 
contributions to the energy appears to be technically harder to obtain than 
for odd p. 

So, for p = 2, 4,... etc and for large meson masses the multiskyrmion 
mass is proportional to the baryon number, and the average (volume) mass 
density decreases as ~ 1/ ^/B. At large B the thickness, or width of the shell 
is given by hh ~ Xj^fm - i.e., it does not depend on the 5-number and the 
mass density in the shell is constant, psheii — contradistinction to 

the case of odd p where it grows with B 

When the meson mass is small, as for pions, we can perform the following 
expansion 

\JoP‘ + 125'(3/b ~ a -f 66'f3/{ba) -I-... (53) 

The main contribution to the mass is then Mq = 2^/a^ at ro = and 

the minimum is reached at ~ 2 X^A^ the massless case jS]. 

Tn this respect there is a direct analogy with the {2 + 1) — D model [laiiaiia, where 
the surface energy density of the rings, representing states of lowest energy, and their 
width, do not depend on the topological number when this number is large. 
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As a result we obtain the following expression for the mass of the multi- 
skyrmion, given the mass term in the Lagrangian, to hrst order in 6': 


M{B) 






m=0 



5'P \ 

2ha‘^ ) 



2 + 


( Qvrtf' 


(54) 


At large values of B the relative contribution of the rn? correction is constant 
(since b' is constant at large B and I/B'^ const, rjB —»• const = 0.3613), 
in contradistinction to the case of odd values of p, and the value of M{B)/B 
from (El is independent of B. Note that the difference of the multiskyrmion 
masses, between the p = 2 and p = 4 cases, is given by: 


M{B,p = 2)- M{B,p = A) 
M{B, m = 0) 




2m?‘ 


i;b(1 - 


(55) 


For the pion mass ~ 0.13, and (El gives the value ~ 0.01, as shown in 
Table H in agreement with the numerical data, for greater m the agreement 
is not so good, since the case of p = 4 is more difficult to describe analytically. 

For any even p the mass term gives a contribution to the multiskyrmion 
mass which is constant at large baryon numbers (relatively), and which de¬ 
creases with increasing p, as then 6' decreases. This is in agreement with the 
numerical results presented in the previous section. 

The radius of the multiskyrmion state, to hrst order in the mass term, 
can also be rewritten as 


ro - ro,^=o(l - - ^o,m=o(l - 


371 , 


(56) 


Since S' decreases with increasing p, the radius of the multiskyrmion in¬ 
creases, in good agreement with the numerical results of the previous section. 
From dH and (El we have also that 


rB{P2) - rsiPi) ^Mb{p2) - Mb{pi) 

_ _ -<_ 

rslPi) ~ Mb{pi) 


(57) 


which is verihed to a good accuracy for B larger than ~ 10. 

To summarize, in the case of even p, he., p = 2, 4..., the multiskyrmions 
have the structure of empty shells; the mass and i?-number densities are con¬ 
centrated in the envelopes of these shells and the energy per unit B decreases 
with increasing B, assymptotically approaching a constant value. 
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m 

o 

o 

o 

o 

(approx) 

0.362 

1.300 

4.131 

0.00989 

0.07602 

0.25053 

0.00987 

0.07615 

0.25089 

0.0100 (SMA) 
0.0528 (SMA) 
0.2149 (LMA) 


Table 6: {M{B,p = 2) — M{B,p = 4:))/M{B,m = 0) for numerical solutions 
(column 2 and 3) and according to SMA (IH^ and LMA ()51I52|) . column 4. 


6.3 Even p; ‘the inclined step approximation’ 

A natural question then arises: to what extent the structure of multiskyrmions 
and their properties depend on the parametrization we have used. Of course, 
we have to satisfy the boundary conditions on the prohle function: /(O) = tt 
and f{oo) = 0 and the function should minimize the value of the mass (j22I)- 
However, the prohle / could have been decreasing according to a law which 
is different from (EH), thus giving us different mass and H-number distribu¬ 
tions. But it is just the property of the Lagrangian (1^ that produces the 
above-mentioned bubble structure as this structure leads to a low value of the 
mass. Another, perhaps the simplest, example of a description that we can 
make is provided by the ‘toy’ model of “the inclined step” type jH]. Such an 
approximation is cruder than “the power step” considered previously. How¬ 
ever, it has the advantage that the calculations can be made for arbitrary p. 
Hence, we shall mention it here and compare its results with what we have 
obtained before. 

Let W be the width of the step, and vq - the radius of the multiskyrmion 
state, dehned by the value of r at which the prohle / = 7r/2. Then we can 
approximate the prohle function by / = 7r/2 — (r — rojvr/hL for Tq — lT/2 < 
r < To + W/2. This approximation describes the usual domain wall energy 
(see, e.g. [Hj) to within an accuracy of ~ 9.5%. 

Next, we write the energy in terms of hh, tq (recall that W ~ ro/6 in 
terms of the previous parametrization) and minimize it with respect to both 
these parameters thus hnding the approximate value of the energy. The case 
of p = 1 was considered previously jS], and since the case of other odd p is 
similar, we restrict our attention here to the case of even p’s. 

Thus for an arbitrary even p we have: 


1 m+Wl2 W- 

wLu,2 + iv 


3 2W^ 

p4^2 


(58) 


The volume term ~ rj] is absent, and since hh -C tq at large B, we retain the 
term ~ rglT on the right hand side of (jSHl) and omit other terms. Then, for 
the classical mass of the multiskyrmion, we have (the second and 4-th order 
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terms were presented in |^): 


M{B,ro,W)c^ — 
on 




p2 


( 59 ) 


The minimization with respect of vq is straightforward and it gives us 


M{B,W) ~ — 
Stt 




m2pp2 7r2\V2 


pz 


+ yj 


+ BI- + W 


(60) 


while = p\/^/[AmyJl + p^n"^/ (2m2lT2^j_ 

We can now consider the two opposite cases. In the case of a large mass, 
when 2m?W‘^ ^ we can expand {w?W‘^jp^ + 7r^/2)^/^ ~ mW/p + 

7 r^p/(4mlT) and obtain 


A/(B, W) ^ T [,v + v^m/p) + T 

This gives us 


(61) 


4m 


V 


(62) 


for IT’”*"' = 7r[(l + v^,^sp/(4m))/(l + mA/3,^B/(2p))]^'^^ where 
At large m the simple formula (El provides an asymptotic (at large B) value 
of the static energy per unit B which, as is easily seen, is in a good agreement 
with our numerical results, i.e., for m ~ 4.13 it gives M/B ~ 1.666 for p = 2 
and M/B ~ 1.408 for p = 4 which agree with the numbers in Table El to 
within (3-4)%. 

On the other hand, the large mass limit cannot be assumed when p is 
large and so we can consider only the small mass limit: 2m?W‘^ -C p^vr^. In 
this case we can consider the m? dependent term in EHD as a perturbation, 
as it was done in jH] and in the section above. Now IT*"®" = n and 


M{B) ^• 

The radius ro is now given by P| 


(63) 


rl ~ (3J/8)^/^ (64) 

and so, tq ~ y/B, and the corrections to tq can be easily found, following the 
steps similar to those of jHI . The difference of masses of the p = 2 and p = 4 
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cases is reproduced well; however, this approximation is too crude at higher 
values of p. 

To conclude, we see that the results obtained within “the inclined step” 
approximation reproduce well the results of the preceeding subsection for 
even p and describe well the transition to higher values of p where the small 
mass limit can be applied. Further rehnements and improvements of this 
analytical discussion are possible, e.g. subasymptotics of the 5-number de¬ 
pendence can be calculated, but we shall not do this here since our results 
are already in a good agreement with the numerical data of section 5, and 
as the asymptotic behaviour of the solitons mass also is well understood, 
(idHimm . etc. 


7 Conclusions 

We have investigated different possible forms of the mass term in the SU{2) 
Skyrme model, concentrating our attention on a class of terms involving a 
parameter p. We have found that the case of even p, but p 7 ^ 0, in the 
parametrization of the mass term, ~ (1 — cos(p/))/p^ is of special interest. 
In this case the contribution to the static mass proportional to the volume 
of multiskyrmion is absent, and the multiskyrmion states, at any value of 
the chiral meson mass m in the Lagrangian, are the spherical bubbles empty 
inside, their energy and baryon number being concentrated on the surface 
of the bubble. The energy per 5-number decreases with increasing 5 and 
approaches a constant value, which is proportional to when m is large. 
Thus as 5 increases the shell conhgurations become more bound. 

For odd values of p and p = 0, there is a volume-like contribution to the 
static mass of multiskyrmions, within the rational map approximation, which 
grows faster than 5. Thus it is responsible for the asymptotic behaviour of 
the multiskyrmion mass ~ and makes the multiskyrmions unstable 

with respect to the decay into skyrmions of smaller values of 5. These 
unbound conhgurations have recently been observed by Battye and Sutcliffe 
0, for p = 1, for relatively low values of 5 even for the pion mass. We have 
also observed them for the other odd values of p: p = 3 and p = 5. When the 
mass m is small, the conhgurations decay is into a5 = 17or5 = 7 shell only 
for very large values of 5. As the mass increases, the decay becomes possible 
for smaller values of 5. When p is even but non null, the conhguration can’t 
decay into smaller shells. Thus if we want our states to be bound states we 
should consider p even, or a combination of terms with p even. 

Note that the fact that some rational map conhgurations do not form 
classically bound states suggests that the states of minimal energy for those 
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values of p and B can be of a different form. They might correspond to 
embedded shells, or not have any shell structure at all. This needs to be 
investigated further but this can be done only by solving the full equations 
of the model. 

The conhgurations we obtained for p > 1 have lower energy than those 
considered within the standard mass term and, for this reason, they have 
a good chance to hnd realization in nature, not only in elementary particle 
and nuclear physics, but also in astrophysics and cosmology. Of course, what 
really happens for physical nuclei is unclear as our results are classical; i. e., to 
compare them with physical nuclei we would have had to compute quantum 
corrections, and this has not yet been done for nonzero modes. 

Investigations of multiskyrmions could be extended also to variants of 
the model where the higher derivative terms (six order, eight order, etc.) 
are included into the effective Lagrangian. The studies performed in [12] for 
the case of the 6-th order term and recently in HE] for some generalizations 
of the Skyrme model including an 8-th order term in the chiral derivatives, 
have shown that topological structures of minimal energy conhgurations are 
the same for these model extensions as in the original variant of the model, 
for values of baryon numbers not too large. Therefore, most probably, the 
observations of bound large B states made in the present paper for certain 
modihcations of the mass term will be conhrmed in such generalizations of 
the model, although this requires detailed studies. 
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